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Abstract 

We conjecture the factorized scattering description for OSP{m/2n)/OSP{m — l/2n) super¬ 
sphere sigma models and OSP{m/2n) Gross Neveu models. The non unitarity of these field 
theories translates into a lack of ‘physical unitarity’ of the S matrices, which are instead 
unitary with respect to the non-positive scalar product inherited from the orthosymplectic 
structure. Nevertheless, we find that formal thermodynamic Bethe ansatz calculations appear 
meaningful, reproduce the correct central charges, and agree with perturbative calculations. 
This paves the way to a more thorough study of these and other models with supergroup 
symmetries using the S matrix approach. 


1 Introduction 


The field theory approach to phase transitions in disordered systems has realized major progress 
over the last few years, thanks to an ever deeper understanding of two dimensional field theories. 
Conformal invariance, combined with elegant reformulations using supersymmetry iii , and 
a greater control of non unitarity issues Q ^, |^, now severely constrains the possible fixed points 
&§• In some simple cases, perturbed conformal field theory, combined with the use of current 
algebra symmetries, has even led to complete solutions |0, 0. Some of the models of interest in 


the context of disordered systems have also appeared independently in string theory |1T, IT], and 
more progress can only be expected from the cross fertilization between these two areas. 

Remarkably, the chief non perturbative method, the integrable approach, has not been pushed 
very far to study these models. This is a priori surprising. For instance, several disordered 
problems involve variants of the OSP{m/2n) Gross Neveu model, which formally appears just as 
integrable as its well known 0{N) counterpart. The standard way of proceeding to study such a 
model would be to determine its S matrix, and then use the thermodynamic Bethe ansatz and 
form-factors to calculate physical properties. This approach was pioneered in the elegant papers 
1^, ^j, and revived in jQ, but so far the subject was only touched upon in our opinion; for 


instance, although the S matrix of the OSP{2/2) Gross Neveu model has been conjectured jlTj, 
no calculation to justify this conjecture has been possible. Super sigma models have also been 
tackled, this time in the context of string theory |jl^, but there again results have only been very 
partial, and the S matrix approach even less developed than for super Gross Neveu models. 

The main reasons for this unsatisfactory situation seem technical. While there has been 
tremendous progress in the understanding of the sine-Gordon model and the 0(3) sigma models 
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- the archetypes of integrable field theories - models based on other Lie algebras are only partially 
understood (see |^, 17| for some recent progress), and the situation becomes even more confusing 


when it comes to super-algebras. One of the main difficulties in understanding these theories is 
physical, and related with a general lack of unitarity - a feature that is natural from the disordered 
condensed matter point of view, but confusing at best from a field theory stand point. Another 
difficulty is simply the complexity of the Bethe ansatz for higher rank algebras, in particular 
super algebras. While these equations can be written sometimes (see the recent recent tour de 
force Ip^ ), finding the pattern of solutions - the generalized string hypothesis - is a daunting task 


even for the trained expert []^|. 

Integrable field theories and lattice models go hand in hand, and the foregoing confusion 
seems to extend to spin chains based on superalgebras. Although the formalism is by now well in 
place to write the integrable Hamiltonians, their continuum limit is not well understood. In the 
case of ordinary algebras for instance, it is known that this continuum limit is a Wess Zumino 
model on the group: whether this is true or not for superalgebras has been a matter of some 
debate [20|. Note that in some cases, the super spin chain is better understood than the field 

] in the spin quantum Hall 


theory: this is the case for instance of the s/(2/l) spin chain of 


effect, whose relation to the traditional (super) Yang Baxter formalism is also not understood at 
the present time. 

Our purpose in this paper is to develop the integrable approach for the case of OSP{m/2n) 
field theories. We will discuss two kinds of models, the supersphere sigma-models, and the Gross 
Neveu models, mostly for algebras OSP{l/2n). In each case, we will conjecture a scattering 
theory, whose striking feature will be the lack of unitarity of the S matrices, as a result of the 
super group symmetry. We will argue that formal thermodynamical calculations do make sense 
nevertheless, and illustrate this point for both types of models. 


2 Algebraic Generalities 


There are two basic integrable models with 0{N) symmetry, the Gross Neveu model and the 
sphere sigma model = 0{N)/0{N — 1). Once their integrability is proven, the scattering 

theory is determined by implementing the action of the symmetry on the space of particles, and 
by requiring factorization. This is not always an obvious task, because of issues of bound states 
and charge fractionalization. For instance, the scattering theory for the 0{2P -|- 1) Gross Neveu 
model was completed only very recently |^. However, the scattering of particles in the defining 
representation has been known for a long time [^], and this is where we would like to start here. 

Scattering matrices with 0{N) symmetry can generally be written in terms of three indepen¬ 
dent tensors 

= aiE + a2P + asl ( 1 ) 

where we have set 


nil 


Pjf 

= 


H Jl 

ri2*2 

iiji 

= 

*1 *2 


corresponding to the graphical representation in figure 1. 

We are interested here in models for which none of the amplitudes vanish. Specifically, for N 
a positive integer, there are generically two known models whose scattering matrix for the vector 
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Figure 1: Graphical representation of the invariant tensors appearing in the S matrix. 
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Figure 2: Pole structure of (T 2 as a function of N. 


representation has the form 


), with none of the fii’s vanishing. They are given by 
2iTT 


fji = - 


0-3 = 


{N -2){m-e) 

2i7r 

0-2 


0-2 


{N - 2)6 


with two possible choices for < 72 : 

r{i-25)r(i + 55) r' 


<72 ( 6 ») = 


-N-2 ^ 2iTT 




N-2 


2i7r 


r(s) r(i-5|;)r(i± 


1 


N-2 


_^ 1 F 

2in ' ^ 


1 


N-2 


+ - 2 - 

^ 2*71 


(3) 


(4) 


The factor uj does not have poles in the physical strip for iV > 0, and the corresponding S matrix 
for > 3 is believed to describe the 0{N)/0{N — 1) sphere {S^~^) sigma model. The factor 
<7^ does not have poles in the physical strip for N < 4. For iV > 4, it describes the scattering 
of vector particles in 0{N) Gross Neveu model. Recall that for = 3,4 the vector particles are 
unstable and disappear from the spectrum, that contains only kinks. Some of these features are 
illustrated for convenience in figure 2. 

Note that at vanishing rapidity, the scattering matrix reduces to S'(0 = 0) = This is in 
agreement with the fundamental particles being bosons in the sigma model , and fermions in the 
Gross Neveu model p5| . 

Our next step is to try to define models for which A" < 1, in particular A^ = 0, or A^ negative. 
A similar question has been tackled by Zamolodchikov [^] under the condition that particles be 
“impenetrable”, that is <7i = 0. The (standard) procedure he used was to study the algebraic 
relations satisfied by the objects E, I for integer N, extend these relations to arbitrary N, and 
find objects (not necessarily N x N matrices) satisfying them. In technical terms, the algebraic 
relations turned out to be the defining ones for the Temperley Lieb algebra , for which plenty 
of representations were known. The most interesting N = 0 case (corresponding to polymers) 
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Figure 3: Graphical representation for the defining relations of the BW algebra. 


could then be studied using the 6-vertex model representation. It could also be studied using 
algebras OSP{2n/2n), or algebras GL{n/n). 

In trying to address the same question for models where ui 0, it is natural to set up the 
problem in algebraic terms again. The objects E, P, I can be understood as providing a particular 
representation of the following Birman Wenzl algebra, defined by generators Ei, Pi, i = 1,... 
and relations 


together with 


and 


PiPi±iPi = Pi±iPiPi±i 

p.2 = 1 

[Pi,Pj] = 0, \i-j\>2, 


EiEi±iEi — Ei 

Ef = NEi 

[Ei,Ej] = 0, |i-j|>2, 


PiEi = EiPi = Ei 

EiPi±iPi = Pi±iPiEi±i = EiEi±i 


(5) 


( 6 ) 


(7) 


These relations can be interpreted graphically as in figure 3; operators E define a sub Temperley 
Lieb algebra [|^]. 

The natural extension of what was done say for polymers would be to look for vertex repre¬ 
sentations of the Birman Wenzl algebra. However, this does not seem possible. The point is that 
the full Birman Wenzl algebra has two parameters, and the representation furnished say by the 
spin one vertex model will have, for instance, that Pi ^ This is a property natural from 

the knot theory framework where this algebra comes from, but disastrous for the construction of 
physical S matrices, where particles cannot “go under” another. Extending the definition of the 
S matrix to arbitrary values of Ai thus seems problematic. 

It is easy nevertheless to extend it to negative integer values of A^. Indeed, the Birman Wenzl 
algebras arise from representation theory of 0{N), and most of the properties of these algebras 
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generalize to the superalgebras OSP(m/2n)[^]. Instead of the vector representation of 0{N), 
take the vector representation of the orthosymplectic algebra, of dimensions (m, 2n). For m ^ 2n, 
the tensor product with itself gives rise to three representations. Taking I as the identity, E as 
(m — 2n) times the projector on the identity representation, and P as the graded permutation 
operator (the extension to the case m = 2n is easy), it can be checked indeed that the relations 
(11))®)® obeyed with N = m — 2n. More explicitly, in the usual case, the matrix elements 
of E are obtained by contracting the ingoing and outgoing indices using the unit matrix. In the 
OSP case, they are obtained similarly by contracting indices using the defining form of the OSP 
algebras 

/ 0 0 \ 

J = 0 0 -4 (8) 

V 0 4 0 / 

In formulas, we set i = i,i = 1,... m, i = n + i,i = m + 1,... ,m + n, i = i. We set x{i) = l,i = 
m + 1,... ,m + n, x{i) = 0 otherwise, so p{i) = x{i) + x{i). One has then 


while the graded permutation operator is of course given by 

pi2*2 _ (xi2 d2 
njl V J 11 jl 


(9) 

( 10 ) 


This realization of Birman Wenzl algebras was first mentioned in the very interesting paper 
It thus follows that the natural ortho-symplectic generalization of the S matrix of the 0{N) Gross 
Neveu model (or sphere sigma model) does provide a solution of the Yang Baxter equation, and 
realizes algebraically the continuation to values of N equal to zero or negative integers. Let us 
now discuss how meaningful this can be physically. 

For this, let us recall some basic features about Yang-Baxter versus graded Yang-Baxter. 
In all cases, the Yang Baxter formalism deals with two related objects that are usually called 
R,R in a general context, S, S in the context of scattering theory, and differ by some (graded) 
permutations. 

In the ordinary case, we reserve the unchecked symbol to the matrix obeying Ri 2 {u — 
v)Ri 3 {u)R 23 {v) = R 23 {v)Ri 3 {u)Ri 2 {u — v), where u,v are spectral parameters. The equiva¬ 
lent of this relation for the superalgebra case is the graded Yang Baxter equation, and it involves 
signs [M|: 


— v)R^i}^-^{u)Rf^^ (^y'j(^_l^pih)p{i2)+p{ki)pik3)+pik2)pik3) _ 
(n — u)(— l)^*^*2)p(*3)+P(h)p(A)3)+P(fcl)p(fc2) 


( 11 ) 


where p{k) = 1,0 is the parity of the k coordinate. These signs occur because, in the graded tensor 
product formalism, R 13 acts on the first and third components, hence giving rise to potential 
minus signs when commuting through the elements of the second component. An ordinary 
(super) R matrix does not solve the graded (ordinary) Yang Baxter equation. However, if R does 
solve the graded Yang-Baxter equation, the object R^j = (—solves then the ordinary 
Yang-Baxter equation, so it is easy to go from one point of view to the other. 

In the ordinary case, one can also consider the object R = PR, P the permutation operator: 
this is what we gave in formula (|I]) for the case N a positive integer. It satisfies a different relation, 
Ri 2 {u — v)R 23 {u)Ri 2 {v) = R 2 d,{v)Ri 2 {u)R 23 {u — v). Observe that this relation now involves only 
neighboring spaces in the tensor product, and thus is insensitive to grading. If R were to solve 
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the graded Yang Baxter equations instead, the same relation would be obeyed by the matrix 
R = PR, where now P is the graded permutation operator. Whether R satisfies the ordinary 
or the graded Yang-Baxter equation, it follows that matrices R do satisfy the same equation. 
Conversely, a solution of R12R23R12 = R23R12R23 can be interpreted as arising from a graded 
or a non graded structure. The graded Yang Baxter equation appears more as an esthetically 
appealing object than a fundamental one. It is especially nice because it admits a classical limit, 
and fits in the general formalism of quantum super groups |3l[| . 

In the context of scattering theories, which are our main interest here, it is convenient to 
define the S matrix through the Fadeev Zamolodchikov algebra |^]. Theories based on su¬ 
pergroups will have a spectrum of particles containing both bosons and fermions. Their cre¬ 
ation and annihilation operators will be denoted and obey for instance zj{0i)zj{92) = 

S^j{6i — 92)zj{62)zj{6i). The consistency of these relations requires that S satisfies 
the graded Yang Baxter equation, or, equivalently, that S satisfies the ordinary Yang Baxter 
equation. Amplitudes of physical processes are then derived in the usual way. An important 
feature is that the monodromy matrix, which describes scattering of a particle through others, is 
built out of S like in the non graded case (the same thing happens for integrable lattice models 

ID- 

Taking therefore our OSP S matrix, and the S matrix that follows from it, S = aiE + a 2 l + 
cr^P, it is natural to ask about the physical meaning of these amplitudes. This reveals some 
surprises. Crossing and unitarity are well implemented in the cases when the particles are bosons 
or fermions. Mixing the two kinds does not seem, a priori, to give rise to any difficulty. For 
instance the relation S{6)S{—6) = S{9)S{—9) = I holds in the graded case with proper choice 
of normalization factors. It will turn out however that in the graded case, the S matrix is, as a 
matrix, not unitary Q. It is thus difficult to interpret our S matrices in terms of a ‘physical’ 
scattering. The most useful way to think of the S matrices will probably be as an object describing 
the monodromy of wave functions, like in imaginary Toda theories |16, Q. Crossing follows then 
from S{i7T — 9) = ai{9)I + a2{9)P + a3{9)E, with an obvious graphical interpretation, and charge 
conjugation being defined through the defining form of the OSP algebra. 

Leaving aside the unitary difficulty, the usual formal procedure thus selects once again the 
factors cT^ as minimal prefactors, with the continued values N = m — 2n. The question is then 
to establish the relation what field theory, if any. 

Obvious candidates are the OSP{m/2n) Gross Neveu model with action (in all this paper, 
normal ordering is left implicit) 


S 



IIL lb 2 

+ 'iPrB^P^r + Y + PrBjr + 9 (V’iV’k + BlIr - ^lBr) 


( 12 ) 


where the V' are Majorana fermions of conformal weight 1/2, and the jSj are bosonic ghosts of 
weight 1/2 as well. Perturbative calculations of the beta function im snggest that this model 
behaves like the continuation of the 0{N) Gross Neveu model to the value N = m — 2n. Similarly, 
the natural generalization of the sphere sigma model is a super sphere sigma model, which can be 
described as the coset OSP{m/2n)/OSP{m — l/2n) . There again, pertnrbative beta functions 
do match. It is therefore natural to expect that the S matrices built on OSP{m/2n) will describe, 
depending on the prefactor cr^, these two models in the appropriate physical regimes. This will 
be discussed in the next section. 

^This is a stronger violation of unitarity than in cases like the Lee-Yang singularity, where = 1 still holds, 
but unphysical signs appear in S matrix residues. For a thorough discussion of unitarity issues see 
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3 The 0SP{l/2) sigma model S matrix 

3.1 The S matrix 

To make things more concrete, let us discuss the case N = —1, and its realization using OSP{l/2). 

Instead of the Gross Neveu model, it will turn out to be easier to study the equivalent of the 

( 2 ) 

sigma model, because of its relation with the O 2 Toda theory and spin chain. 

The solution of the graded Yang Baxter equation relevant here is the well known OSP{l/2) 
one, given by 

1 « j 


R 


OSP{l/2) 


1-2,^ 


2271 *- 


e 


P + —I + 

2i7r ITT — 


:E 


where we have chosen the normalization factor for later purposes, I is the identity. Denote the 
basis vectors in the fundamental representation of OSP{l/2) as b, fi, / 2 . The operator E is given 
by the matrix 

/ 1 -1 1 X 

E=l 1 -1 1 I (14) 

in the subspace spanned by (6, 6), (/i, / 2 ), (/ 2 , /i) in that order, E = 0 otherwise. In that same 
subspace, the graded permutation operator reads 


(15) 



The operators E, P satisfy the defining relations of the Birman Wenzl algebra with N 
The non graded R matrix meanwhile reads 


= - 1 . 


R, 


1 


OSP{l/2) 


1-3^ 


T G 

zzvr ITT — 0 


E 


(16) 


Let us now discuss the issue of unitarity. While R{6)R{—6) = R{9)R{—6) = 1, R, R, and 
R as matrices, are unitary only with respect to an indefinite metric induced by the supergroup 
structure. Explicitly, one has for instance 


nbb fy*bb p/ 1/2 p*/l/2 p/ 2/1 p*/2/l _ 1 

^bb^bb - ^bb ^bb ~ ^bb ^bb - ^ 


(17) 


and in fact R conserves a scalar product that allows for negative norm square states < ff\ff >= 


— 1, all others equal to +1. It is well known indeed |36| that the structure of OSP{l/2) is not 
compatible with a positive scalar product. The mere presence of supergroup symmetry leads 
necessarily to the existence of negative norm-square states, and therefore to unitarity problems. 

The resulting scattering matrix is therefore non unitary, in the usual sense. This is a conse¬ 
quence of the orthosymplectic supergroup symmetry, and originates physically in the non uni¬ 
tarity of the field theory described by the S matrix. This does not prevent one from using the 
S matrix at least to describe the monodromy of the wave functions, as we will do in the section 
devoted to TBA. Similarly, this S matrix could also be used to describe aspects of the finite size 
spectrum [^, ^]. 

An intriguing remark is that, although the matrix R is not unitary, its eigenvalues happen 
to be complex numbers of modulus one (the same hold for R and .R)^, and there are reasons to 
^We thank G. Takacs for suggesting this may be the case. 
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believe that this is true for the eigenvalues of the monodromy matrices involving an arbitrary 
number of particles. This means that the non unitarity situation is not as stringent as say in the 

.(b 


case and that, for instance, the spectrum of the theory in finite size will be real. 

Let us now consider the ‘scattering’ theory that is the continuation of the sphere sigma model 
to = — 1: we take the OSP{l/2) realization, and as a prefactor It then turns out that 
the S matrix is identical to the one of the Toda theory for a particular value of the coupling 
constant! This will allow us to explicitly perform the TEA, and identify the scattering theory 
indeed. While we were carrying out these calculations, we found out two papers where the idea 
has been carried out to some extent already: one by Martins [37|, and one by Sakai and Tsuboi 
|38|. Our approach has little overlap with these papers, and stems from our earlier work on the 
theory instead. 

( 2 ) 

To proceed, we now discuss the Toda theory in more details. 


3.2 A detour through 02 

This theory has action 

1 


( 2 ) 


S = ^ I dxdy + {dy^f + A(2e 

OTT ' 


(18) 


The conformal weight of the first held is Ai = ^, while the one of the second is A 2 = /3^. The 

dimension of A is such that [A]^L“^^2-4/ii _ so [A] = L~'^L^~3 i.e. [A] = and 

the “effective” dimension (i.e. twice the conformal weight) of the perturbation is d = 0^. 

The domain we shall be interested in primarily corresponds to > 1. We will parameterize 


= 2 ^ 


(19) 


so 


h^ = 


hi = ^ = [A] = L The case t = 2 corresponds to hi = i, and the limit t —> 00 to 


2’ 


The massless or massive nature of the theory depends on the sign of A and on the value of 
|39|. For < 1, the theory is massive for A < 0, but for the region we are interested in, A > 0 
is required, and we will restrict to this in the following. 

In the t £ [2, 00 ] domain, the scattering matrix has been first conjectured by Smirnov ||40t| . 
The spectrum does not contain any bound states, and is simply made of solitons with topological 
charges ±1,0 (where the topological charge is defined as q = f dx4>)- The relation between 

the mass of the solitons and the coupling constant reads 1^] 


A'^ = - 


r2(/32/4)r(/?2 


167r3 r2(l -/32/4)r(l -/32) 


ttM 


2 


1 3(2-/32) 


\/3r(i/3) 


^ * 3(2-/3±; J 


( 20 ) 


Near fS'^ = 2, which will turn out to be the point with OSP{l/2) symmetry, setting (3'^ = 2 — e, 
one has A^ oc eM^^. 

The S matrix is proportional to the R matrix of the Izergin Korepin model |^. Although 
this may seem laborious, we will write it explicitly here. Introducing the parameter 




vr/?^ 


32-/32 


( 21 ) 
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and the variables X = e 


-27r6»/5g^ p _ gi7r/2gi7r/3$^ write 


5 = Sr 


1 


\5p5 _ _|_ p-1 _ p 


i? 


with [^, ^] 


R\l 

pO,l 

-^1,0 

K-l 

pOO 

-f^OO 

pOl 

-^01 

p—1,1 
“ 0,0 
pi,o 
“ 1,0 


— A ^+p ^—p 


= Ro'p = K-i = ^-T,0 = Ap^ - A 
= Rhil = Xp - X~^p~^ + p~^ - p 


+ p ^ — p^ +p ^ — p + p^ — p ^ 


R 

R 

R 


0,0 

1,-1 

- 1,1 

1,-1 

1,-1 

1,-1 


- 1,-1 
- 1,-1 

1,0 _ 5 - 1,0 _ 

0,1 “ ^ 

1,-1 _ 

- 1,1 “ 

Xp^ - x-^p-^ 

,Rl};o = A(p® -p)+ p~^ - p~^ 

^°’°i = A(/-l) + l-/ 

Ro~l = X-^{p-^ -p~^) +p^ -p 

^0,0^ = A'Hl -P“^) +p~^-l 

X{p^ — p — p^ +p~^) +p^ — p~^ 
X~^{p~^ — p~^ — p + p~^) +p^ — P~^ 


V ^ +P ^ -p^ 


The normalization factor admits the representation 


So = — exp 


'f 


^^^- 3 ia.e/ 7 r Sinh( 3 a;) cosh(3a;(2^ - 7 r)/ 4 ) 


U) 


It is equal to the amplitude for the scattering processes 11 
In the case ^ ^ oo, one checks that 


sinh(3a;^/27r) cosh(3w/2) 

11 . 


1 


-R 


R. 


osp{l/2) 


( 22 ) 


(23) 


(24) 


(with b ^ 0, /i ,2 


\5p5 _ ^-lp-5 _j_ p-l _ p 
±1) up to an irrelevant gauge transformation. Moreover, it turns out that 


39 

2z7r 


1 - 


36 

2in 


0-n 


(25) 


or So = cr^ — cr^ for = — 1 , confirming the identification of the S matrix in the limit 
t —> oo with the OSP{l/2) “sphere sigma-model” S matrix. 

This coincidence has a simple algebraic origin. Indeed recall H, 0 ] , that the 


.( 2 ) 


Toda 


theory has symmetry Uq{a'^'^), q = . The Dynkin diagram for the algebra turns out to 

be almost identical to the one for the algebra osp(l|2)(^) [^], as represented in figure 4, although 
in the latter case, one of the roots is fermionic, and therefore the basic relations involve an 
anticommutator instead of a commutator. 

It can be hoped that for some particular value of g, the ( 7 -deformation of one algebra gives 
rise to the other, and this is what we shall now demonstrate - namely that there is a mapping 
between Uq{a^'^) and I7(osp(l|2)*^^)), for q = i. This should not come as a surprise, and has 
algebraic roots going back as far as |Q. For recent related works, see [^, ^]. 


9 












osp (1/2) 



(1) 



Figure 4: Dynkin diagrams for and osp{l/2)^^\ 


( 2 ) 

Traditionally, the Cartan matrix of is written as 
relations are 

= 0 

\Hi,Ej\ = CLijEj 

[Hi, Fj] = —QijFj 


-4 


-4 


, and the commutation 


-Hi 


[Ei,F,] = 6, 




Qi-Qi 


n-^i 

Qi = 


(26) 


This means in particular that the generators Eq,Fo,Hq satisfy a [7^4(oi) algebra 

[Hq,Eo] = 8Eq 
[H o,Fo] = -8Fo 


[Eo,Eo] — 


qHo _ q-Ho 

- q-^ 


(27) 


while the generators satisfy a Uq{ai) algebra 

[Hi, El] = 2Ei 
[Hi, El] = -2Fi 


[El, El] = 


qH^ _ q-m 


q-q 


-1 


(28) 


-2 


Commuta- 


The Cartan matrix of osp(l|2)(^) on the other hand reads usually ^ ^ ^ 

tion relations are similar to (p6D, but involve anticommutators instead of commutators for the 
fermionic generators. The generators eo,/o,/io satisfy thus a oi algebra 


[ho,eo] = 4:60 
[ho,fo] = —4/o 

[eoi/o] = ho 

while for the generators , V'l, hi one has 

hi,'ijj\ = 

[hi,'t/;i] = --01 

{ip\,ilii} = hi 


(29) 


(30) 
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Taking q = i for ) makes the subalgebra generated by Eo,Fq,Hq a U{ai) algebra. The 

value q = i for the other deformed oi was already observed in |^| to allow a simple relation with 
a fermionic algebra, a fact also used in mapping Ui{ai) onto a supersymmetric J\f = 1 algebra. 
Here, observe that by setting Ei and '0i = q^^^~^^^‘^Ei/{q + q~^), , one finds, for 

representations where Hi is even (the only ones of interest in our case), that, when q ^ i, 


{V’LV'i} = ^ 

in agreement with the anticommutation relation for U{osp{l\2)) if hi = The rest of the 
relations then are in complete agreement, up to some straightforward changes of normalization. 

We conclude that, restricting to representations with Hi even, the two algebras are isomorphic. 
Since this constraint is satisfied in the case at hand, the osp(l/ 2 )(^) symmetry of the Toda 
theory is thus explained. 


3.3 Thermodynamic Bethe ansatz 


Throughout this paper, we will use the thermodynamic Bethe ansatz to calculate physical prop¬ 
erties of our theory. It is a priori unclear whether the method - which involves maximizing a 
free energy - makes much sense in a theory whose Hamiltonian is not hermitian, but the results 
we obtain seem perfectly meaningful, like in other similar examples. Two additional remarks 
about the TBA are relevant. First, the scattering matrix appearing in the auxiliary monodromy 
problem (diagonalizing the matrix describing the effect of passing a particle through the others) 


is not S but S. This means that, although the S matrices of the and a' 2 '^ differ because of 

the grading, the objects used in the TBA (like the S matrices) are identical, and known results 
about 02 ' Toda theories can be used. Second, one may worry that mixing bosons and fermions 
could give rise to problems in applying the TBA. This is not quite so however. Most TBA’s 
known so far - and the ones we will introduce here will be no exceptions - allow at most one 
particle in a state of a given rapidity. As discussed in Zamolodchikov [48|, this corresponds, in 
the diagonal case, to having S'))(0) = —(—1)^, where E is the fermion number of particle i. In 
our case, we have S^-(O) = For the supersphere sigma model, the particles with bosonic 

internal labels i = 1,... ,m will be bosons, so P-t = (—1)^. For the super Gross-Neveu model, 
the particles with bosonic internal labels are now fermions, so = —(—1)^. In both cases, the 
required result holds. 

The TBA analysis can be performed using the well known strategies. The only difficulty is 

the diagonalization of the monodromy matrix, which involves solving an auxiliary problem based 
( 2 ) 

on the 02 ^ vertex model. String solutions for this model were not known before, but they can 


easily be obtained using our recent results on the case, 
equations have the form 


Setting 7 = ^, 


the Bethe 


n 


sinhi(yi - - iy) 

sinhi(yi -Ua + i^) 


n 


sinhi(yi 

sinhi(yi 


yj - 2i'y) sinh^(pi 
yj + 2^7) sinh^(yi 


Vj + n) 
Vj - il) 


(31) 


where the yi are Bethe roots, and the Ua are spectral parameter heterogeneities (corresponding 
to the rapidities of particles already present in the system). The solutions of these equations in 
the thermodynamic limit are as follows. The y’s can be 1, 2,..., t — 1 strings, or antistrings. In 
addition, it is possible to have a t string centered on an antistring, or to have a complex of the 
form y = yr P ^ + m ■ 
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After the usual manipulations, one ends up with equations for the pseudoenergies, that can 
be represented using a TBA diagram. The ‘left part’ of the diagram corresponds to the following 
equations 


T 


>3(0 - e') *ln(l + {5j 


where we denote (f)p{9) = 2 coshfpe/ 2 ) ^ f * di^) = I-oo ~ ^')9{0')- We use in the following 

the Fourier transform 

/(^) = / (33) 

so if * g) = ‘^'^fg, and We introduce the other kernel if defined by ^ 

In addition, there is a set of equations providing a closure on the right part. 


T 


Together with 


- 0') * In (^1 + - (fsie - e') * In (^1 + - 

~ ^0 * In fl + — 0') * In f 1 + )/'^\ 


i=l 


= -MO - o') * In (1 + + MO 

+ + +ifie 


0 ')*ln(l + e-^o(^')/^) 


and 



MiO - o') * In (^1 + + 2M0 

Y^ i’iO — o') * In ^1 + )/'^^ _|_ 

i=l,2 


0 ')*ln(l + e^o(^')/^) 

0') *ln (^1 + 


(34) 


(35) 


(36) 


Finally, the asymptotic conditions eo(0 —> 00 ) —> mcosh0 must be imposed. This system can be 
conveniently encoded in the diagram of figure 5. The free energy per unit length reads as usual 

/ oo la , . 

—m cosh 0 In (1 + ) (37) 

-00 ^ 

We will consider the more general case of twisted boundary conditions, by adding a phase 
factor in the trace Z = Tr ^ q topological charge. The kinks have therefore a 

fugacity (e^*“/*“^, 1). We concentrate on the central charge, which is expressed in terms of the 
quantities x = in the limits of large and small temperature. At large temperature (UV), 

the Xj go to constants, which solve the following system (here we set A = e*“, which appears in 
the equations due to a renormalization of the spin P9||): 

xo= (l + xi)i/2 (^1 + ^) ' 
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Figure 5; Incidence diagram for the TEA of the (anisotropic) theory. Nodes are associated 
with the pseudoenergies e, and the cross indicates the presence of a massive asymptotic behavior 
for eo- 


Xn= (1 + Xn+l)^/^ (^1 + 

Xt_3 = (l + Xt_4)^/^(l+Xa)^/^(l + AXa)^/^(l + A"lxa)^/^ + ^ (1 + Xfo) 

Xa= (1 +Xt-3)^/^ ^ (l + AXa)“^/^(l + A“^Xa)“^/^(l + Xfe)“^ 

Xb= (1 + Xt-s) (^1 + (1 + AXa)“^(l + A“^Xa)“^(l +Xb)“^ (38) 

and recall that there are three (like the dimension of the fundamental representation) nodes with 
a common value of Xa- The solution of this system is 


Xa = 


Xb = 


(?'+l)a • (j+4:)a 


(t-l)a 


j = 0,..., t - 3 


2t 


(t+l)a 

Sin 

2£ 

sin^ 

(£- 1 ) 0 : 


2£ 

TTTr"— 


(39) 


l+Xj 


What we will in general call the UV contribution to the central charge is ci = 

At small temperature (IR), the Xj go similarly to constants solving the same system but with one 
less node on the left, because xq ^ 0. We first consider the case a = 0, i.e. periodic boundary 
conditions for the bosons, antiperiodic boundary conditions for the fermions. In that case, the 
UV sum of dilogarithms gives a contribution (t — 1), while one gets a similar contribution from 
the IR sum after t ^ t — 1, C 2 = t — 2. The central charge is thus c = ci — C 2 = 1, as expected. 

[Here we include two specialized remarks: 

1 ( 2 ) 

A point of some interest is 7 = corresponding to /ii = ^. In that case, the ' Bethe 
equations do coincide (after a shift y —>■ y + y) the Bethe equations that appear in 

g2 

solving the sine-Gordon model with This point is in the attractive regime, with one 

soliton and one antisoliton of mass m, and one breather of the same mass. It is easy to check 
that in that case, the TBA is in fact identical with the well known SG TEA indeed. The 
equivalence between the two theories is not so obvious when one looks at the actions. 
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Also, as is related to a 2 ^\ so does the a' 2 '^ theory bear some resemblance to the Toda 


theory with the following action 


5 = — 
Stt 


dxdy 


+ {dy4>if - A' ^ 

.i=l ^ . 


(40) 


Here the perturbation has a single dimension h = j3‘^, and the dimension of the coupling is 


[A] = Parameterizing (3 in (^0[) by 0^ = turns out that the free energy of the 

02 ^^ theory is exactly half the free energy of the theory, once the fundamental masses have 
been matched. This fact does not appear obvious in the least when one compares perturbative 
expansions!] 

( 2 ) 

Twistings and truncations of the model are of the highest interest and have been widely 
discussed in the literature 42, 44|. Twisting (that is, putting a charge at infinity) in such a way 


that becomes a screening operator of weight A = 1, gives the central charge c = 1 — 

RSOS restriction is then possible for t even, giving rise to the minimal model The 

perturbation in the minimal model has then weight A 21 = 1 ~ ^ (its coupling is real, and the 
sign does not matter because it has only even non vanishing correlators). Meanwhile, the lowest 
weight Ai 2 = becomes negative for t > 4, after which the effective central charge reads 

Cg£f = c — 24 /ii2 = 1 — t{t\) • twist in such a way that becomes a screening 

operator, giving the central charge c = 1 — RSOS restriction is then possible for t odd, 

giving rise to (•j_i )/2 perturbed by the operator of weight A 15 = 1 — f • 0 

The twisting can also be studied with the TEA using now a 0 0. The UV sum of dilogarithms 
gives then a contribution (t — 1) — 3^^^, while one gets a similar contribution from the IR sum 
after t —> t — 1. The central charge is thus c = 1 — 

In general, twisting terms affecting nodes ‘far to the right’ of the TEA diagram do not affect 
the central charge in the isotropic limit. Indeed, if a were to remain finite here as t ^ 00 , the 
central charge of the twisted theory would be still c = 1. We shall however be interested in giving 
antiperiodic boundary conditions to the kinks of charge q = ±1, which translates into a phase 
that blows up like a ~ tvr as t ^ 00 . As a result, the central charge of interest is c ^ — 2 , in 
agreement with the sigma model interpretation to be discussed next. 


Finally, we notice that choosing a = 27r leads to xt _4 = 0, and a truncation of the diagram 
to the one represented in figure 6 . This is the same as folding the TEA for the RSOS model 
with central charge c = 2 — The first model in the series has = |, the next one 

Ceff = I (the latter TEA has some fascinating properties, due to the fact that | + | = 1). This 
was first observed in [ |50| . We will comment about the relation of these models to OSP{\/2) in 
the conclusion. For TEAs related with in other regimes, see [^, 


3.4 The OS'P(l/2) limit, and the relation with the sigma model 

As explained previously, the OSP{l/2) scattering theory can be studied by taking the t ^ 00 
f2) 

limit of the model. The identification could in fact be seen directly by identifying Eethe 

^Notice that the combination 2 — x = |,resp. | for t even (resp. odd). In fact, the perturbative series for the 
free energy always has the same structure, and does not exhibit parity effects as t is changed. But the physical 
interpretation does, and rightly so, since for (^ 21 , only even correlation functions do not vanish, while for (^ 15 , all 
correlation functions are a priori non vanishing. 
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1 


t-6 


t-5 


f 2') 

Figure 6 : Incidence diagram of the truncated TEA. This TEA also describes perturbations 
of the OS'-P(l/ 2 )( 4 _t)/ 5 'C/( 2 )p_ 4)/2 models for t even, see later. 


( 2 ) 

equations. This seems a bit strange at first, because the equations do not have a structure 
that is reminiscent of the osp{l/2) Dynkin diagram. One has to remember however that the 
osp{l/2) Eethe ansatz equations are peculiar, and their structure is not related with the Cartan 
matrix in the usual way. They read in fact [^] 


n K — fj^a — i _ TT ~ ~ -r-r Aj — Xj + i 

Aj — Poi T ^ Aj — A^ + 2i Aj — Xj — i 


( 2 ) 

and match the equations in the t —> oo limit, with y = yA, u = y/i, 7 —> 0. 
The Toda theory (|i^) can then be rewritten in terms of a Dirac fermion as 


S = 


(fix 

'2^ 


+ '^\pi>L + A 




(41) 


(42) 


the perturbation is the sum of a term of dimension hi = ^, and a term of dimension /12 = 2 
(the relative normalization between the two fermionic terms is irrelevant, since it can be adjusted 
by XipR, —> A^^V'rj oi' similarly for left fermions). It is likely that this model could 

be directly diagonalized using the coordinate Eethe ansatz, like the ordinary massive Thirring 
model, but we have not carried out such a calculation. Conserved quantities can be found in 
terms of the fermions; the first ones are {dips'll: r)... 

The twisted theory meanwhile has c = —2, Cgg- = 1, and the perturbations both acquire 
dimensions (1,1). This can be identified with a symplectic fermion theory with action [^] 

-S' = y ^ (43) 


Here 771 and 772 are two fermionic fields with propagator, in the free theory, < 771(2;, .2)772(0) >= 
— hizz. Notice how non unitarity is manifest in (^) as well as (^). 

From the point of view of the twisted theory, the perturbation involves two fields of weights 
(1,1) which should be identified with (f) 2 i and 4)i^ respectively, using fusion relations. That 
both fields appear is not unexpected, since the c = —2 point is a limit, and should have the 
characteristics of both t even and t odd. 

The identification of 9^771772 with (p 2 i can actually be completed accurately, by comparing 
the four point functions as calculated in the fermion theory and the minimal model using the 
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Dotsenko Fateev general results [M|. An interesting sign subtlety appears in that case. Indeed, 

di_ir]idi_ir ]2 = dr]idr ]2 + dr]idr] 2 , and if we call this operator O, < 0(l)0{2) >= —because 

^ 12^12 

of anticommutation relations. Hence, should actually be identified with i4>2i- In fact, 

( 2 ) 

when one compares the amplitude of the perturbation in the Toda theory and the twisted 


version |39|, one finds that, with the usual normalizations, A positive gives rise to the amplitude 
of (j) 2 i being purely imaginary, that is the coefficient A' in real. The sign of A' is irrelevant, 
as only terms even in A' will appear in the perturbative expansions of physical quantities. 

It would be very interesting to complete the identification of </>i 5 with but we 

have not finished this calculation. Note however that there is little doubt this identification is 
correct, as there is no other object with the right dimension and statistics in the symplectic 
fermion theory. Defining O = riirj 2 d^r]id^r] 2 , one finds < 0(1)0(2) >= ) . The massive 


perturbation with is obtained with a coefficient that is real and positive near 0^ = 2 |39]. 
Therefore, we expect A" in (|4^ ) to be positive. Note that the apparition of logarithms in the two 
point function of the perturbing operator makes the field theory (^) a bit problematic. Issues of 
renormalizability arise in particular, and it is probably better to think of (^) as a sector of (^ ) 
rather than the defining theory. This is reflected in the structure of the TBA: although c = —2 
can formally be obtained as the UV value of the central charge in the untwisted model, this value 
appears only after proper analytic continuation of the dilogarithms. Indeed, the fugacity given 
to end nodes of the TBA diagram is e*“, and as as t ^ oo, a ~ tn, it winds an infinite number 
of times around the origin : in practice, following the free energy would presumably require 
following analytic continuations on an infinity of different branches, a difficult task at best. 

The fermions can always be rescaled to bring the action into the form 


S = 


cPx 




(44) 


where again the coupling A is positive. We will now see how this related to the super sigma 
model. 

In general, the coset space OSP{m/2n)/OSP{m — l/2n) has dimensions (m — l,2n) and 
can be interpreted as the supersphere g£ interest here is m = n = 1, and 

corresponds to the S^’'^ supersphere, parameterized by the coordinates 

1 

xi = 1 - -r]i'n2 

6 = vi 
6 = 112 

such that Xi + ^ 1.^2 = 1- The action of the sigma model will generally be of the form 


(45) 


1 / ™ 

- ( 


+ 


(our convention is that the Boltzmann weight is e~^). The beta function will be to first order 
/3 oc — (m — 2n — 2)g'^, so for the region m — 2n < 2 in which we are interested, the model will 
be free in the UV and massive in the IR for a negative coupling constant, g = —\g\. In the S^’'^ 
case, this action therefore reads 


S = - 


1 


Px 


df,gidf,g2 - 


(46) 
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Note that a rescaling combined with a relabeling can always bring this action into the form 


S = J d^x 
,( 2 ) 


d^riid^ri2 + 


matching the t —> oo limit of the theory, with A oc |( 7 |. 


(47) 


4 Supersphere sigma models and integrable superspin chains 

The relation we uncovered between and OSP{l/2) extends immediately to the case of 
and OSP{l/2n): one can establish, for general values of n, the relation between the quantum 
affine algebras, the Bethe ansatz equations, the scattering matrices etc. We thus propose that 
the S matrix with OSP{l/2n) symmetry, represented in (|^),(^),(|^ with N = 1 — 2n, and the 
prefactor , provides an analytic continuation of the 0{N)/0{N — 1) “sphere” sigma model to 
this value of N. 

Of course, the analytic continuation of the sigma model should be interpreted as the coset 
OSP{l/2n)/OSP{0/2n). The effective central charge of the UV limit is Cgg- = n, while its true 
central charge will be c = —2n. For the ordinary sigma models, the UV central charge is — 1, 
so the UV value in the analytic continuation just matches. 

The Toda theory has an interaction term of the form 

g^/3(01-</-2) g^/3(</>2-03) ^iV2l3<l>n 

The dimension of vertex operators exp(i ^ dj(j)j) is h = ^ — (5o X) (5j, where ho measures the 

twist, and the central charge is c = n(l — 12hQ). The operators in the interaction term all have 

h = 1 when ho = and P = y/2. In that case, c = —2n, while all the operators and 

dimensions ( 1 , 1 ). 

The manifold relevant for OSP{l,2n) on the other hand is i.e. a purely “fermionic 

sphere”. For instance, can be parameterized by 

1 / ,1 

3^1 = 1-2 idiV2 + mm) - -^mmmm 
6 = hi 

6 = m 

= m 

= m (48) 

The action of the sigma model is not particularly illuminating; it involves four and six fermions 

couplings, and reduces to 2n symplectic fermions in the UV limit. Like in the n = 1 case, it can 

(2) 

be matched onto the appropriate limit of the theory. 

( 2 ) 

On the other hand, it is also possible to extend the analysis of the ' TBA to arbitrary value 
of n, so we also know the TBA for this scattering theory, which is simply given by a Z2 folding 
of the TBA. The TBA is represented in figure 7. 

Notice that there are n massive particles: while for N integer positive the 0{N)/0{N — 1) S 
matrix has no bound states, with simply N fundamental particles (in the vector representation), 
poles do enter the physical strip for N < 2. For the value N = 1 — 2n we are interested in, the 
masses of the particles are rrii oc sin 2 ^^, i = 1,..., n. The UV central charge is easily checked 
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Figure 7: Incidence diagram for the TEA of OSP{\/2n) sigma models (the diagram has n rows). 


to be Cgg- = n. We do not know how to obtain the central charge of the untwisted theory, as this 
would require a knowledge of the ‘closure’ of the TEA diagram for twisted theories, an unsolved 
problem when n > 1. 

Our results have an immediate application to the study of quantum spin chains. Indeed, the 
Eethe equations which appear in the solution of the OSP{\/2n) sigma models are similar to 
the ones appearing in the solution of the integrable OSP{l/2n) chains studied in particular by 
Martins and Nienhuis |^. More detailed calculations show that these chains are critical, and 
that they coincide at large distance with the weakly coupled supersphere sigma models, that is, 
a system of 2n free symplectic fermions. This is disagreement with the conjecture in |20, 37| that 
this continuum limit should be a WZW model on the supergroup: although the central charge 
agrees with both proposals, detailed calculations of the thermodynamics or finite size spectra 
show that the WZW proposal is not correct, and confirm the sigma model proposal instead. A 
similar conclusion holds for OSP{m/2n) when m — 2n < 2. That the spin chain flows to the 
weakly coupled sigma model is certainly related with the change of sign of the beta function when 
m — 2n crosses the value 2, but we lack a detailed understanding, similar to the ones proposed 
, 56], of the mechanisms involved. 


m 


5 The super Gross Neveu models 

5.1 Generalities 

If we consider a scattering matrix defined again by i) ,(^), but now with the prefactor a 2 in¬ 
stead, it is natural to expect that it describes OSP{m/2n) Gross Neveu models, the analytic 
continuation of the 0{N) GN models to 0{m — 2n). Having a control on the diagonalization of 
OSP{l/2n) scattering matrices will allow us to study this scattering theory easily, and confirm 
the identification for these algebras. Notice that since the 0{N) scattering matrix has no poles 
in the region N < 2, the roles of the GN and sigma models are completely exchanged in the 
domain of values of N we are considering. 
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The OSP{m/2n) Gross Neveu models read 



/3i97i + + 9 

i=l j=l 



(49) 


This theory has central charge c = ^ — n, effective central charge Cgg- = ^ + 2n. The beta 
function for this model is of the form I3g oc (m — 2n — 2)g^, the same as the one for the 0{m — 2n) 
GN model. For m — 2n > 2, it is thus positive, so a positive coupling g is marginally relevant 
- this is the usual massive GN model - while a negative one is marginally irrelevant. If instead 
we consider the case m — 2n < 2, these results are switched: it is a negative coupling that 
is marginally relevant, and makes the theory massive in the IR |^. The case m = 1 should be 
described by the foregoing scattering theory. 

Note that the GN model is equivalent to the appropriate WZW model with a current current 
perturbation. Indeed, the system of m Majorana fermions and n symplectic bosons constitutes 
in fact a certain representation of the OSP{m/2n) current algebra. The level depends on the 
choice of normalization; it would be called A: = 1/2 in 


k = —1/2 in |58], k = 1 elsewhere. 
We adopt the latter convention here, and thus the level k WZW model based on OSP{m/2n) 
has central charge 

(m — 2n)(m — 2n — l)k 

(50) 


c = 


2{m — 2n — 2 + k) 


Particular cases are OS'P(0/2n), which coincides with the SP{2n) WZW model at level —k/2, 

and OSP{m/Qi), which coincides with the 0{m) WZW model at level k. Super symmetric 

space theorems give rise to free helds representations at level A: = 1, where c = at level 

; n n 1 (m—2n)(m—2n—l) sdim OSP(m/2n) , . ii, 4 . j. j.- 

k = m — 2n — 2, where c = - - = - 2 “-^^—-—INotice that the representation 

at level —2 for OSP{2/2) described recently in is a particular case of the super symmetric 
space theorem discussed in Goddard et al. (for A: = — 1 in their notations). 

The OSP{m/2n) Gross Neveu models present additional non unitarity problems not encoun¬ 
tered in the sigma models discussed above. To tackle these problems, we first discuss the simplest 
case of all. 


5.2 The OS'P(0/2) case. 

We consider the case of the GN model for N = —2, corresponding formally to O5P(0/2), i.e. 
a (3^ system. The S matrix should act on a doublet of particles, and reads, from the general 
formulas 


„ , TT iO 

5 = tan ( - — 


r (4 - 1 - r (_ 

-^ 12 ^ 2in) V 2*71/ 


/ -1 
0 


'1 _ e_\ 

V 2 2i7r ) 


r(B7) 


0 

1 

6—in 


V 


0 - 
0 


e 

0—in 

0 


e 

6—in 
1 


0 

0 

0 


6—in 

0 -1 / 


(51) 


It turns out that S = f tanh (| — ^) SsciPsG ~ where Ssg is the soliton S matrix of the 
sine-Gordon model. At coupling /3gQ = Svr, it coincides with the S matrix of the SU{2) invariant 
Thirring model, or the level 1 WZW model with a current current perturbation. 

'*In 1^, the four fermion coupling is defined through combinations -1-7/)-'/’+ = + '/’£'/’«) i so what 

is called g there is the opposite of our convention. 
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Figure 8; Incidence diagram for the TEA describing the current-current perturbation of the 
(anisotropic) /37 system (or O5P(0/2)) model. 

The scattering matrix is thus the same as the one for the k = 1 SU{2) WZW model up to a 
CDD factor. This CDD factor does not introduce any additional physical pole, but affects the 
TEA in an essential way. 

Note that the N = —2 Gross Neveu model can also be considered as a current current 
perturbation of the SU{2) WZW model at level k = —or the SU{1, 1) model at level k = ^. 

To study the TEA, it is useful as in the sigma model case to consider the anisotropic defor- 
mation, with the sine-Gordon part now corresponding to = -j-. The diagram is represented 
in figure 8. The UV solutions have to obey the usual SG equations plus the fact that 

xo = + + (52) 

The solution is obtained by setting Xj = {j -|- a)^ — 1, j = 0,..., t — 3, xt -2 = xt-i = t — 3 + a, 
and letting a oo. The contribution to the central charge in the UV is then ci = t, the number 
of nodes. The solution in the IR is obtained by discarding the first node, and then coincides with 
the usual IR solution of the SG equations, with a = 1. The contribution to the central charge is 
equal to C 2 = t — 2, the number of nodes minus one. The final central charge is thus c = 2, as 
expected for the effective central charge of the j3^ system. 

The same calculation with a fugacity e *-i gives c = 2 — This is because in the UV, 

all the x’s are still infinite, giving rise this time to ci = t while in the IR, the re’s are the 

same as the ones for the ordinary Sine-Gordon model, with C 2 = t — 2 — 6|5f This result 
requires explanations; in particular, setting a = (t — l)7r and letting t ^ oo as in the sigma model 
case gives c = — oo! 

5.3 The role of zero modes 

We want to consider in more details the (5^ system with action 

•S’ = y ^ {hd^/L + PRd'yn) (53) 

The propagators are jr{z)Pr{w) = —PR{z)'yR{w) = We can ‘bosonize’ the ghosts by 

introducing a scalar field ^ = (j)R + (pR, such that (Pr{z)(Pr{w) = — \n{z — w). We also introduce 
fermionic ghosts r]R{z)CR{w) = ^Riz)r]R{w) = and thus 

IR = 'yi = e~'^^r]L 

Pr = e~'^^d^R, Pl = (54) 
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The corresponding action is then 


S = — [ cPx = TT f (9<h5<h (55) 

Hit J 2tt J 

The P'y Hamiltonian is 

^ J dx {fSidx'jL + PrOxIr) (56) 

with commutators [/3 l(x), 7 ^( 1 /)] = ^(^(x - y), [/?ij(x), 7 ij(y)] = -^(5(a: - y). The 17(1) current 
is given by Jr = 'JrPr = dcpR, Jr = 'JlPl = -dcpL- The topological charge is 

Q = '^ J dx {Jr - Jl) = ^ J dx dx^ (57) 


The topological charge of 71 , and 'Jr is 1, while the charge of Pr and Pr is —1. 

A key feature of this system is the existence of zero modes. With periodic boundary conditions, 
it is indeed easy to see that [H, f JxPr^r] = [H, f Jxjr^r] = 0. It follows from this that, if we 
add to the Hamiltonian a term of the form —hQ, the system will fill up with an infinity of zero 
mode particles of /3 or 7 type depending on the sign of h, sending the ground state energy to — 00 . 
The theory is thus unstable without a mass term. In the current algebra language, the infinite 
dimensional space associated with the zero mode decomposes into lowest weight representations 
of SU{1, 1 ) 1/2 of ‘angular momentum’ j = —| and j = —The conformal weight of these states 
is A = — |, giving the effective central charge Cgg = 2 for a c = — 1 theory indeed. 

The mass term (which is actually a current current perturbation) in the OSP{0/2) GN model 
does stabilize the theory. To see how, let us add to the action a term 

' 55 ’ = -^ J (fx{-iRPR - 'JrPr) + ^ y d^x {'JrPr - 'JrPr^ ( 58 ) 

The classical minima of S + 6S occur for 'yR = jr = c and Pr = —Pr = b, and, turning to the 
Hamiltonian formalism, the minimum energy becomes then 


If 

L ~ 27t g 


We now recall the RG equation for the coupling constant g in (^): 


9 = - 25 ' 


(59) 


(60) 


From this, the coupling constant at scale 1/m goes like g = 2 in Lt/m ' constant term is a UV 

cut-off, provided here by the field h. It follows that 

yEgs = --h^ ln{h/m) (61) 

L TT 

at leading order. If m —> 0 (// —> 0), we recover the result Egg —00 anticipated before. 

We will comment more on the behavior of the OSP{0/2) and other GN models later. For the 
moment, our goal is to explain the behavior of the central charge in the anisotropic case obtained 
in the previous section. So, we now consider the case where an anisotropy is imposed on the 
system by adding a coupling of the form Jr Jr. More explicitly, consider 

d A = J (fx{jRPR - 'JrPr) - j d^x 2-fRPR'yRPR (62) 
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It is easy to calculate the ground state energy at leading order as <7 —> 0, which turns out to be 
finite now, even though the theory is still massless: Egs/L = — Anisotropy has stabilized the 
UV theory. 

We can now calculate this ground state energy using the S matrix approach. To do so, we 
perturb the action by a term of the form In the bosonized version, this reads 

e~^'^drjLVLd^CRdCR + e^’^drjRTjnd'^^LdCL 

The anisotropic term changes the kinetic term to ^(1 — We can renormalize the field 

so the kinetic term looks as before, and then the exponentials in the perturbation become 
with = T^. Non local conserved currents are then obtained using 

exp (^ 1 ^) d'^^Rd^R, exp d’^riRrjR 

of dimension Ac = 3 — 2//3^. They lead to a quantum deformation of the algebra with 
quantum parameter q = Setting q = e *-i, this corresponds to a thermodynamic Bethe 

ansatz diagram with t nodes (including the source one), i. e. the TBA studied in the previous 
section and represented in figure 8 . The point is, we now have the correspondence between the 
coupling g in the anisotropic action and the parameter t in the anisotropic TBA, with, at small 
g or large t, g ^ j. To use this TBA, we finally need to establish the correspondence between 
the magnetic field and the kinks fugacity: setting where T is the temperature, 

gives h = Using this, the TBA result c = 2 — f does match the ground state energy 

Egs = — at leading order as s' —> 0. We thus have explained the TBA results of the previous 
section in the light of the ground state instability of the f3'y system. 

For the [3j system itself, we of course obtain c = — 00 for any non trivial fugacity of the kinks. 
This can also be understood as follows: the partition function of the /Sy conformal field theory 
with periodic boundary conditions is infinite because of the presence of a bosonic zero mode j^ . 
On the other hand, in the periodic sector, the TBA gives a hnite result, with the central charge 
c = 1 + 2. Therefore, the TBA approach must describe a renormalized partition function, divided 
by the infinite contribution of the bosonic zero mode. In the antiperiodic sector where there is no 
such zero mode, the result of this division is to give zero, or, formally, a central charge equal to 
— 00 . Nevertheless, the dependence of the ground state energy on a can be predicted and checked 
using the TBA, which provides another non trivial check of the S matrix. 

5.4 The OSF(l/2n) case 

We now get back to the OSP(l/2n) case. The TBA turns out to have a simple description in terms 
of 02 ^ again. Consider therefore, not the SU{2n + 1) GN model, but a related scattering theory 
with only two multiplets of particles, corresponding respectively to the defining representation 
and its conjugate. Considering more generally the case of SU{P) models, the relation between 
the SU{P) GN scattering theory and this new theory is similar to the relation between the 0{P) 
GN model and the 0{P)/0{P — I) sigma model 13. We will thus call this scattering theory 
‘sigma model like’, but we are not aware of any physical interpretation for it. The TBA equations 
can be written following the usual procedure. They read 

P—1 00 

b=l 1=1 
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Figure 9: Incidence diagram for the TBA describing the anisotropic version of the OSP{l/2) GN 
model. 


(-al 

T 


p-i 


Ia^^4>P * In ^1 + — 0p * In + e 

6=1 

- (<5ai + Sa,P-i)4>P * In (^1 + 


for the pseudoparticles, and 

£io 

T 


ep-1,0 


m cosh 6 

f 

m cosh 6 


p-i 


+ ^ </>ia * In [1 + 

a=l 
P-1 

+ ^ (/)p-i,a * In 


a=l 


(63) 


(64) 


In these equations again, (f)p{6) = </>p-i,a = <?^>i,P-a, and (^ia(a;) = 

This TBA is in fact quite similar to the one of the N = 2 supersymmetric SU (P) Toda theory 
(the generalization of the supersymmetric sine-Gordon model for 517(2)): the difference 
affects only which nodes correspond to massive particles, and which ones to pseudo particles. 
As a result, the central charge is easily determined, c = 2P — 1. Getting back to the particular 
case P = 2n + 1, we can then fold this system to obtain (see the appendix for the proof) 
the TBA for the OSP{l/2n) Gross Neveu model, whose effective central charge reads therefore 
Cgff = i( 2 ( 2 n + l)-l) = 2 n + |. 

As an example, we can discuss in more details the case of the OSP{l/2) GN model whose 
TBA is represented in the figure 9. We will even consider an anisotropic generalization of this 
TBA, where the diagram is truncated to the right in a way that is equivalent to what happened 
in the sigma model case. 

There is now a single color index, and we relabel €aj = ej. Introducing Xj = the 

( 2 ) 

equations in the UV are the same as for the ^ Toda theory, except for the first one that reads 
now simply 

xo = 1 + xi (65) 

( 2 ) 

The closure equations in particular are the same as for the anisotropic model. 

The solution in the case of a vanishing twist first is xj = ^ in the limit a —> oo, 

for j = 0,..., t — 3. In addition one has Xa = and Xb = As a —> oo, all the x's 
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go to infinity but Xa which goes to one. As a result, the UV contribution to the central charge 

isci=t — l + 3x^ = t + ^. In the IR, the modification due to the cj 2 factor is not seen any 

( 2 ) 

longer, and the x's obey the same equations as for the ' case, with a contribution t — 2 to the 
central charge. It follows that c = |, as expected. 

The twisted TEA follows from similar principles as in the sigma model case. This time 
however, because of a = oo, the UV values xj are unaffected by the twisting. The contribution 
to the central charge is c = t — l + ^ + {L\ + L^-i) {xa = 1) = t + \ — The IR values do depend 
on a, with formulas identical to the sigma model case, and a contribution c = t — 2 — 3|5f 
The resulting central charge is c = | — 3^^^^^;^. The dependence on a is similar to what we 
observed in the O5P(0,2) case, for similar reasons. The factor 3 in this formula, as opposed to 
the factor 6 in the O5P(0, 2 ) case, has its origin in the different relations between the physical 
anisotropy and the parameter t in the TEA: in the OSP{l/2) case, 5 ~ 


6 Finite field calculations 


To give further evidence for our S matrices, we now present some finite field calculations. The 
idea, which has been worked out in great details in other cases [ 6 ^ ], is to compare S matrix 
and perturbative calculations for the ground state energy of the theory in the presence of an 
external field. The S matrix calculations are considerably simpler than the TEA ones because, 
for a proper choice of charge coupling to the external field the ground state fills up with only one 
type of particles, with diagonal scattering, and the Wiener-Hopf method can be used to solve the 
integral equations analytically. 

The S matrix calculations are very close to the ones already performed for the 0{N) sigma 
model and the 0{N) Gross Neveu models. In fact, in the region m — 2n > 2, the calculations 
are identical, since the S matrix elements are obtained by continuation N = m — 2n, and, in the 
domain m — 2n > 2, the integral representations are obtained by the same continuation as well. 
For these cases, one thus immediately checks that the continuation of the S matrix to N = m — 2n 
matches the beta functions of the sigma or Gross Neveu models, which are, too, obtained by this 
continuation. 

Things are more interesting in the case m — 2n < 0 in particular, to which we turn now. 
We consider first the OSP sigma model, with S matrices determined by a~^. If we couple the 
external field to a charge of the form 

Q oc J{xidtX 2 - X 2 dtxi)dx ( 66 ) 

the ground state fills up with bosonic particles of the form |1 >+i |2 >, with diagonal scattering 
s = crj + (T^ {9). If meanwhile we couple the external field to a charge of the form 

Q (X J{^idt^2 + ^2dt^i)dx (67) 


the ground state fills up with fermionic particles of the form |I >+|2 >, with diagonal scattering 

s = 

Let us consider this latter case. Using formulas given in the first section, one finds 

+ + _ r(i + x) r(i/2 - x) r(i/2 + a + x) r(i + a - x) 

2 3 r(i-x)r(i/2 + x)r(i/2 + A-x)r(i + A + x) ^ ’ 


24 




where x = -^ and A = j^z 2 - This turns out to coincide with 


______ r(i + x) r(i /2 - x) r(i /2 - a + x) r(i - a - x) 

“ r(i - x) r(i/2 + x) r(i/2 - a - x) r(i - a + x) 

after a continuation A ^ — A. Similarly, 

___ _ r(i + x) r(i/2 - x) r(i/2 - a + x) r(-A - x) 
“ r(i - x) r(i/2 + x) r(i/2 - a - x) r(-A + x) 


does coincide with 


^ r(i + x) r(i/2 — x) r(i/2 + A + x) r(A — x) 

= r(i - x) r(i/2 + x) r(i/2 + a - x) r(A + x) 


(69) 


(70) 


(71) 


after the same continuation. 

This means that, in a TEA calculation, the ground state energy of the OSP{m/2n) sigma 
model coupled to a fermionic charge follows from known calculations about the 0{N) Gross 
Neveu model after formally setting N = m — 2n and performing a continuation iV — 2 —> 2 — A. 
From expressions in we find therefore, for m — 2n < 0 


E{h) - F(0) 


h? 11 /I \^lnln(/i/m) 

27r ^ A - 2 ln(/i/m) \A-2y \T?(h/m) 

1 Cn / lnln(Vm) Y 

N—2\T?(h/m) )_ 


(72) 


with Cn = InF ~ ~ Iv^) 2 + 1, and N = m — 2n. From this we deduce Q the 

ratio of the first two coefficients of the beta function as ^ 


Similarly, from TEA calculations, the ground state energy of the OSP{m/2n) Gross Neveu 
model follows from known calculations about the 0{N) sphere sigma model after formally setting 
N = m — 2n and performing a continuation A — 2 ^ 2 — A. From expressions in |^] we find 
therefore for m — 2n <2 


E{h)-E{Q) = {N-2)^ 


ln(/i/m) — 


1 


A-2 


lnl„(Vm) + OA. + 0(lAA/”> 

\ m{h/i 


im] 


(73) 


where Dn = — ■jy^ln2 — + 77 ^) ~ InF Trom this we deduce the ratio ^ = 

— jy^ 2 - Observe that the leading term follows from the calculations for the (3^ system described 
in the previous section; all that has to be changed is the beta function for the coupling g in 
(]5^) , resulting in Egs = \n{h/m) indeed. Remarkably, it is the ground state of the Gross 

Neveu model that has a leading /i^ln(/i/m) dependence, while the ground state of the sigma 
model has a leading pure h? dependence: the roles of Gross Neveu and sigma model are therefore 
switched compared to the usual 0(A) situation. 


The calculation with a coupling to the first kind of charge (66) for m — 2n < 2 or the second 
type of charge ( |67|) for m — 2n>2 poses difficulties, as the kernel does not factorize in the usual 
way then, so the Wiener Hopf method does not seem applicable. 
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Now recall that for the usual 0{N) sphere sigma model, ^ and for the usual 0{N) 


Gross Neveu model 


N ^ m — 2n, as desired]^. 


’ ~ W^' ratios we found are thns the analytic continuations to 


7 Conclusions and speculations 

To conclude, although more verifications ought to be carried out to complete our identifications, 
we believe we have determined the scattering matrices for the massive regimes of the OSP{m/2n) 
GN and the OSP{m/2n)/OSP{m — l/2n) sigma models in the simple case m = 1, based on 
algebraic considerations as well as thermodynamic Bethe ansatz calculations. 

It is tempting to expect that at least some of our results generalize to other cases OSP(m/2n) 
for m > 1 and m — 2n < 2. In all these cases, we expect that the S matrix of the sphere sigma 
model will be obtained from the conjecture at the beginning of this paper, with N = m — 2n, for 
N < 2. The S matrix of the GN model is probably more complicated. Recall that in the case 

> 2, it is given by the general conjectnre only for N > A. When N <2, we think it is probably 
given by the conjecture only for iV < 0 . 

Observe now that for the nsual 0{N) case, the factors and do not exhibit poles and 
are equal for N = 3,4. For these values, the (unique) S matrix based on the general conjecture 
(|^ describes correctly the sigma model. As for the Gross Neveu model, its description is more 
subtle; it turns out that the vector particles are actually unstable, and that the spectrum is made 
of kinks only. 

In the case < 2 of interest here, the factors and similarly do not exhibit poles and 
are equal for A" = 1,0. These cases would correspond for instance to OSP{3/2) and OSP{2/2) 
respectively. It is very likely that there again, the S matrices describe the sigma model, and not 
the Gross Neveu model, for which the proper particle content has still to be identified. 

The OSP{2/2) case is particularly intriguing. The 5^’^ sphere can be parameterized by 

xi = cos (j){l - ^riir] 2 ) 

X2 = sm4>{l - ^r]ir]2) 

6 = Vi 

6 = 112 (74) 

The action of the sigma model now reads 

5* = -|^ / (Px [(5/,4>)^(l - r]ir] 2 ) + d^r]id^r ]2 - mil 2 dtiVidt^V 2 ] (75) 

where $ is compactified, <I> = $ mod 27r. A rescaling and a relabeling brings it into the form 

S = j cPx [-((9^$)^(1 + \g\'nir]2) + df,r]idf,r ]2 + \g\viV2diimd,,m] (76) 

with now <I> = <hmod -^S=. We see in particular that in the limit (7 —> 0, the action is simply the 
Visl 

one of a free uncompactified boson and a free fermion, of total central charge c = — 1 , and that the 

®The existence of different sectors in the OSP{m/2n) models does not spoil this conclusion, as the beta functions 
are only trivially affected by the twists. 
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boson has a negative coupling : this system therefore coincides with the standard “bosonization” 
of the /dj system in the limit \g\ —> 0 ! 

Notice that as soon as m > 1, the negative sign of the bosons coupling for OSP{m/2n) sigma 
models will have to be handled carefully; presumably the OSP{2/2), jSj example will provide a 
good example of how to do this. 

Finally, we discuss the boundary value N = 2, which exhibits some exceptional features. 
Indeed, from the integrable point of view, the solution of the Yang Baxter equation (combined 
with crossing and unitarity) based on the generic S matrix is not unique for = 2 (in contrast 
with the other values of Y) but admits one continuous parameter 7 . This solution is close to the 
sine-Gordon solution, and is related to it by 


S = a 2 +<73 

St = CTl + 72 
Sr = 7i + 73 

where S, St, Sr are the usual sine Gordon amplitudes 


(77) 


St = -i 


.sinh(87r0/7) 


5 = 


— 


1 


sin(87r2/7) 
.sinh(87r(z7r — 0)/j) 


sin(87r2/7) 


Sr 


Sr = -sm{87r^/j)U{e) 


IT 


(78) 


and U{0) is given e.g. in [^|. 

In the case of 0(2), the existence of this parameter corresponds to the fact that the 0(2)/0(l) 
sigma model, or the 0(2) GN model are actually massless critical theories, the coupling g being 
exactly marginal. The S matrices then provide a massless description of these theories. Since 
72 = 72 , the 0(2)/0(l) sigma model and the 0(2) GN model coincide; their identity follows 
from bosonization of the massless Thirring model into the Gaussian model. The free parameter 
in the S matrices is related with the coupling constant in either version of the model. (Note that 
the S matrices can also be used to describe some massive perturbations. These, however, give rise 
to different type of models than the ones we are interested in, like the massive Thirring model. ) 

It seems very likely that similar things occur for OSP(2n + 2/2n) models as well. The identity 
of the sigma model and the GN model in that case is not obvious, but one can at least check 
using our general formulas that the central charge and the effective central charge do match, 
Ceff = 3n + 1, c = 1. There are on the other hand strong arguments showing that the beta 
function is exactly zero [^, so these models should have a line of fixed points indeed [^, in 
agreement with the S matrix prediction. 


Besides completing the identifications we have sketched here, the most pressing questions that 
come to mind are: what are the S matrices of the Gross Neveu models for non generic values 
of N, what are the S matrices for the multiflavour GN models, what are the S matrices for the 
orthosymplectic Principal Chiral Models? We hope to report some answers to these questions 
soon. As a final related remark, recall that there is an embedding ^ 

OSP{l/2)_2k ^ SU(2), X (79) 

®The level —2k in this formula stems from our conventions; it would be k if it were defined with respect to the 
sub SU{2). 
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and that the branching functions of the latter part dehne a Virasoro minimal model, with 


2 k 


Cosp — 


Csu2 — 


2k+ 3 
3k 

k + 2 


Cvirasoro — 16 


{k + lf 


{k + 2){2k + 3) 


(80) 


For k an integer, the situation is especially interesting. The Virasoro models which appear 
there have p = 2k + 3,q = k + 2-, they are non unitary, and their effective central charge is 
Cgg = 1 — (fc+ 2 )( 2 fc+ 3 ) • These models can thus be considered as UOSP/SU coset models. Their 

perturbation by the operator (^21 with dimension h = 1— coincides with the RSOS models 

( 2 ) 

defined in section 3 as truncations of the a\ theories with t = 2A: + 4. We thus see that the 


supersphere sigma model appears as the limit A; ^ 00 of a series of coset models |W], just like 
the ordinary sphere sigma model say appears as the limit of a series of parafermion theories, 
this time of SU{2)/U{1) type. There are many other interesting aspects of OSP coset models in 
relation with the present paper which we also plan to discuss elsewhere. 
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A Folding the SU{P) TBAs 


The quantization equations for the SU{P) GN model have been written for instance in |]A]: 

P—1 00 

27rPa0 = rua coshp + "^ -kpaj, a = 1 ,...,P-1 

b=i j=i 

P—1 00 

2TTPaj = * Phi 


b=l 1=1 


(81) 


Using Fourier transform; f{uj) = f{(5), one then has 


y(P) 

^ab 


6ab - e'^l 


sinh((P — a)Lo) sinh(6w) 
sinh(Pa;) sinh(a;) 


a >b 


(82) 


with Yab = Tfea- The kernels are logarithmic derivatives of the scattering matrix between 
top components in each of the fundamental representations. The group structure is encoded in 
the densities p of (massless) pseudoparticles, which appear in the solution of the auxiliary Bethe 
system diagonalizing the monodromy matrix. 
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We now set P = 2n + 1. The folded equations then read 

n oo 

27rP„o = ma cosh /3 + ^ * Pbo - X] Paj,a = I,... ,n 

6=1 J=1 

n oo 

6=1 l=l 

For instance, = —e(~’^+^/^)l^l check that the corre¬ 

sponding kernel coincides with a'^ — for N = 1 — 2n, and with the corresponding S matrix 

element in the scattering theory |^. The couplings between pseudoparticles can also be 

( 2 ) ( 2 ) 

checked to arise from the structure of solutions of the a^n Bethe equations, generalizing the 
case. 

The ‘sigma model like’ equations for SU{P) are based on a hypothetical scattering theory 
with physical particles in the vector representation and its conjugate only. They read 


27rPio = m cosh /? + z[P -k pio + * pp_i,o - ^ -k py 


(P) 


doo) 


i=i 


27rPp_i,o = m cosh P + Z[ j_^k pio + Z},]^ p_-^k pp_ifik pp_ij 

i=i 

P—1 OO 

2TTpaj = (<5al + 6a,P-l) * PaO “ 

6=1 l=l 


(84) 


where z[^'^ = Z^p\p i = e and Z^l , = Z^p\. = e 1“^! kernel 

11 P—1,P—1 sinh(Pa;) I 5 P—1 P—1»1 sinh(Pa;) 

is the logarithmic derivative of what is called in 0’ = ZYP. Setting 

again P = 2n + \ and folding gives now 


2pPio = m cosh P + [Zp + * pi,o - ^ o'j * Pij 

i=i 


2TTPaj = * PaO - + ^apnll-h) * hi 


(00) 


6=1 l=l 


-(2?i+l) jp{2npl) 


(85) 


The kernel + z[^YY — ^^iinh^^n+i)^^ ^ differs from the previous kernel + 

Yjkh~^^ by — ^ , which coincides with the Fourier transform of the ratio — for N = 1 —2n. 

l:2ri. -2 cosh cr^ 

Some integral representations to finish (used in the domain < 0) 

6 - 2 “^ - 1 cLo 


lncj 2 = 


10 


^i{2-N)Puj/iT _|_ ^-{2-N)uj^-{2-N)iuiP/n 


g-{2-N)uj _|_ I 


UJ 


and 


i„T= r 

a 2 J-i 


J{2-N)Puj/-it 


cosh doj 

cosh ^ 


( 86 ) 

(87) 


I sinh 0 —2 sin_ 

where ^ 


sinh 0+2 sin 
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